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We examine the thermopower Q of a mesoscopic normal-metal (N) wire in 
contact to superconducting (S) segments and show that even with electron- 
hole symmetry, Q may become finite due to the presence of supercurrents. 
Moreover, we show how the dominant part of Q can be directly related to the 
equilibrium supercurrents in the structure. We also discuss the thermopower 
arising due to an anomalous kinetic coefficient which is finite in the presence 
of supercurrent and in some situations gives the dominant contribution. In 
general, a finite thermopower appears both between the N reservoirs and the 
superconductors, and between the N reservoirs themselves. The latter, how- 
ever, strongly depends on the geometrical symmetry of the structure. The 
paper includes a detailed analytical derivation of the results and an exact nu- 
merical solution of the quasiclassical equations in a few sample geometries. 
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1. INTRODUCTION 

Prom time to time in the last century, the connection between the ther- 
moelectric effects and superconductivity has attracted interest. 1 As early as 
1927, Meissner found that the thermopower is absent in a steady-state super- 
conductor. This is due to the fact that the thermoelectric current is in this 
case counterbalanced by the supercurrent. However, in 1944, Ginzburg 2 
pointed out that there can be other types of thermoelectric effects in su- 
perconductors. Interestingly, many of them are still unobserved 1 , at least 
partially due to the fact that they are based on the same electron-hole asym- 
metry that produces the thermoelectric effects in nonsuperconducting sys- 
tems 3 , where Mott's law connects the thermopower to the tiny difference of 
the effective masses between electrons and holes (i.e., to the nonlinearity of 
the quasiparticle dispersion relation above and below the Fermi sea). Such 
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effects are governed by the factor ksT/Ep, and are thus very small in the 
low temperatures required for superconductivity. Even after six decades, the 
understanding of the experiments on some of these effects is still lacking. 4 ' 5 
For example, the measured thermoelectric flux in a bimetallic superconduct- 
ing ring was orders of magnitude larger than predicted by the conventional 
theories. 

Another effect of thermoelectric type, finite even in the presence of com- 
plete electron-hole symmetry, was observed and successfully explained in the 
turn of 1980's: 6 ~ 9 it was found that the combination of a temperature gradi- 
ent and supercurrent could make rise to a charge (branch) imbalance, i.e., a 
difference in the electrochemical potentials of the normal and superconduct- 
ing components. The origin of this effect is similar to that considered in this 
paper, but it manifests itself in quite a different form in the systems studied 
here. 

In the 1990's, many groups started to investigate the superconducting 
effects induced in normal-metal (N) wires from nearby superconductors (for 
review, see Refs. 10-12). For example, if the N metal is connected to two 
superconducting segments, one can drive a supercurrent through it . It 
was found that this superconducting proximity effect influences both the 
electrical 14 ' 15 and thermal 16 conductances of the N wire, in general breaking 
the Wiedemann-Franz relation between the two. 17 In most cases, these effects 
account for some tens of percent changes in the conductances. 

In proportion, a much greater effect was experimentally observed in 
the thermopower 18-22 of "Andreev interferometers", structures where the 
normal-metal wire was connected to two ends of a broken superconducting 
loop. The observed thermopower was found to be orders of magnitude larger 
with the proximity effect (below the superconducting critical temperature 
Tc) than without it (T > Tq). A thermally induced voltage was found both 
between the ends of the normal metal 18 ~ 20 > 22 ("N-N thermopower" Qnn) as 
also between one end of the normal metal and the superconducting contacts 
("N-S thermopower" Qns)- 21,22 Moreover, it was found that this induced 
thermopower oscillates with the flux applied through the superconducting 
loop. 

It was shown in Ref. 23, including one of the present authors, that Mott's 
law is not in general valid in the presence of the proximity effect, indicating 
that the reason underlying the observations is not necessarily dependent on 
electron-hole symmetry. Indeed, the experimental observations were partially 
explained, without invoking electron-hole asymmetry, in the linear regime for 
Qns by Seviour and Volkov 24 and another effect was suggested for the pres- 
ence of Qnn by Kogan, Pavlovskii and Volkov. 25 The full nonlinear regime, 
including important effects coming from asymmetries in the structure, and a 
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phenomenological explanation that connects the thermopower to the temper- 
ature dependence of the supercurrent, was described by the present authors 
in Ref. 26. This theory seems to be in quantitative and qualitative agreement 
with most of the observations. The aim of the present paper is to explain 
this theory in detail and describe a related effect, thermopower due to an 
anomalous kinetic coefficient, which is the main source of thermopower in 
certain geometries. 

Further, Ref. 27 shows that the presence of supercurrent can lead also to 
a Peltier- type effect. This effect has not yet been experimentally observed. 

This paper is organized as follows. In Sect. we detail the quasi- 
classical formalism applied in the calculation of the thermopower. This for- 
malism assumes complete electron-hole symmetry, and the predicted effects 
are hence not limited by the factor kBT/Ep- The analytic approximations 
necessary for obtaining an expression for the thermopower are presented in 
Sect.lHU including symmetry arguments for the full (without approximations) 
solutions, a phenomenological picture of the effect, and a calculation of the 
thermopower in the quasiequilibrium limit. The general behavior of this 
thermopower is detailed in Sect. ^ and effects related to the geometry of 
the studied structures are discussed in Sect. OD Most of the results in these 
sections are obtained numerically without any approximations (beyond the 
quasiclassical diffusive limit). Section |f3 includes a comparison of the the- 
ory to the experiments and to the other theories on the effect. The main 
spectral coefficients responsible for the thermopower, the spectral supercur- 
rent and the anomalous coefficient Tj are described in Sect. E] Finally, the 
results are shortly discussed in Sect. |^ An impatient reader should capture 
mostly the kinetic equations Q in Sect.EJ skip the detailed analytic deriva- 
tions in Subs. COD and proceed with the main results in Subs. IXn"l and their 
characteristics in Sees. l4.|fT1 



2. FORMALISM 

Many well-understood phenomena in inhomogeneous superconductivity 
have been explained through the use of quasiclassical Green's functions. This 
is in particular the case in mesoscopic superconductivity (see the reviews in 
Refs. 10,12). The essence of the quasiclassical theory is ignoring the short- 
scale oscillations in the relative coordinate, and concentrating only on the 
center-of-mass coordinate of the space-dependent Green's functions. Such 
an approximation ignores the quantum-mechanical interference effects arising 
from short scales of the order of the Fermi wavelength (structural variations 
on this scale can be included through boundary conditions, see below and 
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Ref. 28), and the effects related to the nonlinear parts of the quasiparticle 
dispersion relation (the theory is based on the linearisation of the dispersion 
relation around the Fermi energy). Therefore, it cannot describe the normal- 
metal thermoelectric effects which arise from the energy dependence of the 
effective mass. 

In this paper, we concentrate on the diffusive limit, where all the impor- 
tant length scales are much greater than the mean free path. This is typically 
the case in metallic wires. To describe nonequilibrium effects (finite voltage) 
in the nonlinear regime, we adopt the Keldysh formalism. 29 In the considered 
limit, the Keldysh-Green's function is described by the Usadel equation. 30 



2.1. Keldysh Usadel equations 

The Keldysh-Usadel equation for the Green's functions G(f, E) in the 
Keldysh <8) Nambu space can be written in the static case as 

V • {DGVG) = [-iE 1 ® f 3 + A + S in , G] . (1) 

Here, D is the diffusion constant of the metal, S; n contains the inelastic 
scattering self- energies, and A is the pair potential matrix. The pair potential 
is self-consistently related to the solutions of Eq. JU, 12 but it vanishes in 
normal metals, on which we concentrate here (superconductivity is induced 
by a boundary condition from bulk superconductors). Moreover, G has the 
matrix structure 

G = (°* . G A = -f 3 (G*)tf3 , G K = G R h - hG A , (2) 

where G R , G A and G , the retarded, advanced and Keldysh components, 
are 2x2 matrices, and the last relation is due to the normalization G 2 = 1. 
The electron-hole distribution function matrix h introduced here has two 
free parameters and it can be chosen of the form h = /l1 + frfs, where fx 
and /l are the symmetric and antisymmetric parts, 

f T (E) = 1 - f(fi s -E)- f(fi s + E), f L (E) = f(/i 5 -E)- f(jM S + E) (3) 

of the electron distribution function f (r, E) and the energy E is given with 
respect to the potential /is of the superconductors. 

The Green's function obtained from Eq. (QJ may be used for evaluating 
the observable energy and charge current densities jq and j c : 

3Q = ^ / E 3L dE > jc= Tej h dE ' (4a) 

j L ee i-nr[( Tl I) GVG] , j T ee iTr[(n f 3 ) GVG] . (4b) 
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Here, jl and jt are the spectral current densities, and a is the normal-state 
conductivity of the metal. 

It is useful to parameterize G so that it automatically satisfies the nor- 
malization. Since G 2 = 1 implies that (G R ) 2 = 1, one can parametrize G R 
in terms of two complex functions 9 and x as 

* R ( cosh(0) e^ sin h(0)\ 

" V-e~ ix sinh(#) - cosh(#) J ' 1 ' 

and hence G" 4 and G K depend only on 9, x, /t and To find the equations 
that determine 9 and one ca n consider the retarded part in Eq. (£Q), and, 
neglecting inelastic effects (S; n = 0), obtain the spectral Keldysh-Usadel 
equations in normal metals (A = 0) 

DV 2 9 = -2iE sinh(0) + ^D(V X ) 2 sinh(2#) , (6a) 
DV ■ JE = 0, JE = ~ smh 2 (9)V X ■ (6b) 

These equations describe the behavior of the superconducting proximity ef- 
fect: the pairing amplitude can be written as F = sinh(#) exp(zx), which 
shows that 9 is related to the strength of the proximity effect and x to the 
phase of the superconducting order parameter. 

Using the same assumptions as above and taking the traces Tr[(ri®l) • ] 
and Tr[(ri®T3) • ] of Eq. one obtains the kinetic equations for the spectral 
energy and charge current densities 

DV-j L = 0, jL = V L Vf L -7Vf T + jsf T , (7a) 

DV-j T = 0, 3T = V T Vf T + TVf L +jsf L . (7b) 

These describe the kinetics of electrons and holes in the presence of the 
proximity effect, and, moreover, imply that ji and jx are conserved. The 
spectral coefficients in Eqs. (JJJ are determined from the solutions to Eqs. ||JjJ: 

V L = -Tr[l-G R G A ] = i(l + | cosh^| 2 - |sinh0| 2 cosh(2Im[x])), (8a) 

V T = -Tr[l -G R t 3 G A f 3 ] = -(1 + | cosh#| 2 + | sinh#| 2 cosh(2Im[x])), (8b) 

T= -Ti[G A G R r 3 ] = -| sinhfl| 2 sinh(2Im[x])), (8c) 

js = ^Tr[(G R VG R - G A VG A )f 3 ] = Im[- sinh 2 (0)V X ] = Im\j E ]. (8d) 

Here, T>l and T>t are the spectral energy and charge diffusion coefficients, 
and js is the spectral density of the supercurrent-carrying states. 31 The 
cross-term Tis usually small but not completely negligible. 



P. Virtanen and T. T. Heikkila 



To simplify the problem, one often assumes that wires are quasi- ID 
structures, i.e., much longer than wide, and translationally invariant in their 
latitudinal directions. Thus, one needs only to consider longitudinal variation 
in the quantities, and the gradients may be replaced with ID-derivatives. 



2.2. Boundary conditions 

The Keldysh-Usadel equations cannot handle changes in the structure 
which occur in a distance short compared to the Fermi wavelength or the 
mean free path, so boundary conditions describing interfaces and junctions 
need to be derived using different techniques. 

To obtain the boundary conditions at nodes of one-dimensional wires 
with clean metallic contacts, one may use the the matrix circuit theory. 
28 It yields a Kirchhoff-like conservation law for the spectral matrix cur- 
rents AaGd x G, where A is the cross-sectional area of a wire. In the 9- 
parameterized form, this results in 

^2 AmVOi = 0, A ^Xi = 0, (9a) 

i i 

A i a dL =0, A i a ^T = °> ( 9b ) 



and to the continuity of 9, %, /l and fx- In Eqs. ((HI, the sums go over the 
wires ending in the node. 

The properties of connections between wires and terminals are also de- 
scribed by a matrix conservation condition. 28 In the limit of a clean metallic 
contact, at the interface to a normal reservoir, Green's function G and thus 
also 6>, X: II and fa get their bulk values: 

9 = 0, x = arbitrary, (10a) 




(10b) 



where T is the temperature of the reservoir and /j, its potential, compared 
to fj,s- Thus, near normal metal reservoirs or far from superconductors, the 
coefficients ((HI) obtain their normal-state bulk values T>l — ► 1, T>t — ► 1 and 
T — ► 0. 

At a clean metallic contact to a superconducting reservoir, bulk values 
are again obtained, with the exception of /l for energies inside the super- 
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conducting energy gap A: 

A 

9 = artanh— , x = <t>, h = 0, (11a) 
\E\ > A: /L = tanh(^), (lib) 
\E\<A: j L = 0. (11c) 

Here, <j> is the phase of the superconducting order parameter. We assumed 
here no charge imbalance, fj,s = 0, hence fx = 0. Moreover, the boundary 
conditions for 6 and x imply that near a superconductor interface T>i — > 
and T — ► for energies \E\ < A, which in fact implies that jl = 0. 
Thus no contribution to the energy current density jq through the interface 
comes from inside the energy gap — this is a known property of the Andreev 
reflection. 

With these boundary conditions, the spectral equations © depend only 
on the phase differences and the geometry of the system, while kinetic equa- 
tions depend also on the temperatures and the potentials of the reservoirs. 
The system is, however, fully decoupled from the temperatures of the super- 
conductors for the energies \E\ < A. 

Note that in wires connected to a normal reservoir, the spectral su- 
percurrent je = 0, due to the boundary condition 8 = 0. Since generally 
sinh(#) in these wires, the conservation of spectral supercurrent js im- 
plies the conservation of the phase, thus x = const (x). Similarly, in wires 
connected to a superconducting reservoir, there is the boundary condition 
jl = for \E\ < A, and one can eliminate /l from the kinetic equations. 



2.3. Energy and length scales 

Scaling the spatial dimensions in the spectral equations (JHJ by the length 
L (where L can be one of the wire lengths, or most naturally for the supercur- 
rent, the distance between the superconductors), one finds a natural energy 
scale: the Thouless energy 

f 13 M V e 1 D/(200 
— L 2 \ 0.15 K k B J (L/fim) 2 ' [ ' 

This energy is the inverse time of flight through a diffusive wire of length L. 

In practice, there are two different types of energy scales describing the 
behavior of the spectral functions: Et and the superconducting energy gap 
A of the electrodes connected to the normal-metal wire. 31 The previous is 
important in the case of long junctions L»(o (i.e., Et <C A) and the latter 
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for L < £ {Et It A). Here £o = \/^D /2A is the superconducting coherence 
length. For example, for T <C Tq, in Al £o ~ 200 nm, and in Nb £o ~ 30 nm. 

In the long-junction limit A 3> Et, the energy scale for the solution 
to (|6|7|) with the boundary conditions (|9I11|) depends only on Et- Thus, by 
changing Et, the results may be scaled to fit all systems with similar ratios 
of wire lengths, cross sectional areas and conductivities. 



3. THERMOPOWER IN ANDREEV INTERFEROMETERS 

A temperature difference AT between two points in a metal may in- 
duce a charge current I c . Disconnecting the system from its surroundings 
makes the current vanish, and a potential difference AV arises instead. The 
thermopower Q is the ratio of this potential difference to the temperature 
difference, Q = (AV/AT) I . In this section, we show how the supercurrent 
flowing in Andreev interferometers connects the temperatures and potentials, 
inducing a finite thermopower. 



3.1. Thermopower without superconductivity 

An estimate for the thermopower in the absence of superconductivity 
can be found using the semiclassical Boltzmann equation as, for example, 
in Ref. 3. Such a description can be obtained from the microscopic theory 
presented above, provided that the length scales important for the problem 
are much greater than the Fermi wavelength. 29 This analysis leads to the 
Mott relation for the thermopower 



ir 2 k%T dlna(E) 



dE 



(13) 

E — Ef 



Here, a(E) describes a generalized conductivity — ct(Ef) gives the actual 
conductivity in the linear-response regime — and the logarithmic derivative 
depends on the energy dependence of the elastic relaxation time and the exact 
structure of the Fermi surface. The previous can typically be neglected, and 
the latter is due to the small asymmetry between the dispersion relations for 
electrons and holes. Both of these effects are neglected in the quasiclassical 
theory. There are also other major contributions to the thermopower, e.g. 
the phonon drag, but for simple metals, they should vanish at high and very 
low temperatures. 32 

We can use the temperature dependence in Eq. (|T3*|) for an extrapola- 
tion from the room-temperature values of Q to its value at subkelvin tem- 
peratures. Usual magnitudes for the thermopower at the room temperature 
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Figure 1. The system under examination. The terminals S are supercon- 
ducting reservoirs, N normal-metal reservoirs and they are connected by 
normal-metal wires with lengths Lf., cross sections and normal-state con- 
ductivities cjfc, k = 1 . . . 5. Temperatures T and potentials fi = eV in the 
reservoirs are as shown. Moreover, all contacts are assumed clean metallic. 
Positive directions for currents are indicated with arrows. 

are of the order for pure metals. To obtain a rough estimate at a low 
temperature, we extrapolate 10 from room temperature to a tempera- 
ture of 0.5 K (typical for the studies on the proximity effect), which yields 
Q ~ 10 S^L. Measurements in Ref. 21 gave a stricter estimate Q < 1 ^ for 
their case. 

However, the Keldysh-Usadel equations predict the presence of a large 
thermopower Q ~ 4*-, even under complete electron-hole symmetry, as su- 
percurrents couple the temperatures to the potentials (through the coeffi- 
cients js and Tin the kinetic equations (0). 



3.2. Modeled structure 

Here we focus to study thermopower in the structure shown in Fig. ^ 
where two superconducting (S) reservoirs with a phase difference <f> are con- 
nected through a normal-metal link to two normal metal (N) reservoirs. 
When measuring thermopower in this structure, no current flows in wires 
1 and 2, 

3d = 0, 3d = 0, (14) 

and definite voltages V\ and Vi are induced in the N reservoirs. Here we may 
define several thermopower-like quantities: 

V 2 - Vi Vi V 2 

QNN = 7f 7f ) <?NS,1 = 7f 7f- i QNS,2 = 7f 7jT ■ (15) 

±2 — J- 1 ±2 — J-l 1-2 — 1-1 

Clearly, their behavior with varying parameters need not be similar. That is, 
generally the potential induced in the colder electrode may differ from that 
induced to the hotter electrode. 
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The thermopower appearing in the structure may be estimated by ap- 
plying the formalism described in Sect. E] We have solved Eqs. I)6|7|) numer- 
ically, and derived analytic approximations for the solutions of the kinetic 
equations 0. In the following, we first discuss the analytic results, as they 
provide some insight to the problem, and then examine interesting features 
numerically (Sees. El and EJ . 



3.3. Qualitative picture 

The effect of supercurrents on the thermopower can be understood phe- 
nomenologically as follows. If T\ ^ T2, the temperature-dependent 13 equi- 
librium supercurrent Ig(Ti) in wire 3 is different from Is{T2) in wire 4. (For 
this qualitative picture, we approximate these wires to be at the tempera- 
tures Ti, T 2 of the normal reservoirs.) Conservation of currents should still 
be maintained, so a compensating effect must arise. Should the normal reser- 
voirs be kept at the same potential as the superconductors, a quasiparticle 
current I qp oc (Is{T\) — Jg^)) from them to the superconductors would bal- 
ance the difference. However, when no current is allowed to flow in wires 1 
and 2, a compensating N-S potential difference Vn — Vs oc R(Is(Ti)—Ig(T2)) 
is induced instead. 

The qualitative picture gives the following predictions: First, the in- 
duced potentials should oscillate with the phase difference <p between the S 
electrodes, similarly to the supercurrent. Moreover, the potentials induced 
in the N electrodes need not be the same, especially if the resistances of the 
wires are not symmetric at the left and at the right. Thus, the supercurrent- 
effect also contributes, as in Eq. (|31b|) . to the thermopower Qnn between 
the N electrodes. Such a contribution should be present even in structurally 
symmetric setups, as the resistances are temperature-dependent due to the 
proximity effect. 14 



3.4. Phase behavior 

If the phase difference between the superconducting elements vanishes, 
the quasiclassical equations predict a vanishing thermopower Q (as shown 
below). In the presence of a phase difference, the supercurrent and the 
term T couple the two types of distribution functions, inducing a finite Q. 
Although the solution to Eq. Q cannot be found in closed form, we can 
obtain information about the exact phase behavior of the thermopower using 
symmetry considerations. Suppose we have the solutions x> fh and /t for 
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the Keldysh-Usadel equations l|6l7p in the system. Then, define 

X = -X, = 6, f L = fL, f~T = -f T . (16) 

First, we see that these new variables are solutions to Eqs. (0- This implies 
that T>t = T>t, T>l =T>l,T= — Tand js = —js> so the new variables are also 
solutions to Eqs. Q. They also satisfy the spectral current conservation JjjJ, 
with jt = —Jt and jl = jl, and the boundary conditions (jlOpilj) . but with 
(ft = —(ft and /2 = —fi. Thus, there is a solution with inverted charge currents 
I c , potentials \x and phases (ft. For the thermopower this implies that 



= -Q(-cft). (17) 

lc=0 

Moreover, Q is necessarily 27r-periodic in phase, so this implies that it should 
also be antisymmetric with respect to (ft = nir for all n £ Z. 

The result also implies that in a system with no phase differences or no 
superconducting parts, the thermopower vanishes. This is in agreement with 
the Mott-law prediction due to the quasiclassical approximation. However, 
note that the reasoning above assumes that the condition lfTl)l determines the 
solution uniquely and that the thermopower is well defined. This may not 
be the case in some pathological geometries in the presence of the proximity 
effect, as the condition lfH|l is determined only for the energy-integrated 
currents. In the cases studied in the present manuscript, the convergence of 
numerics hints that the solutions are unique. 



Q{4>) 



AV 



AT 



AV 



AT 



3.5. Analytic approximations 

Although we cannot solve analytically the spectral equations in the 
structure of Fig. ^ we can obtain useful approximations starting from the 
kinetic equations. 

As a first step, we assume that the induced potentials are small compared 
to the temperature, \fj,\ <C fegT. Since fx is related to the potentials, we can 
reason that /t is small, and hence neglect the terms —Td x fT and jsfr in 
Eq. (|7ajl . Comparison with numerics shows that this is a good approximation. 

The kinetic equations, where inelastic effects are neglected, can now be 
integrated over a wire of length L. This relates the currents to the distribu- 
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tion functions fi, fr at the ends of the wire: 

■ MB ~ /l '°> Ma 1 

3L = L(V?) ■ (18a) 

jT =MB^m +jsfm (18b) 



LHVJ 1 )^ 1 ) \Jo V L V T J J V L (x')V T (x 



Here (T>l)t) — T fo ^lJt ^ x are ^ ne dimensionless spectral energy and 
charge resistances. All the coefficients can be calculated numerically by 
solving the spectral equations. We see that Eq. (|T8|l is linear in the dis- 
tribution functions, so we can use the continuity of fx and fi at the nodes 
and write the Kirchhoff-like spectral current conservation equations l|9b|) for 
the structure in Fig. ^ as: 

f 5 (o)\ 



aiAiX.! [ ) + ^3X3 I ^0 > + <tbA 5 X 5 I f ^ ) = , (19a) 
a 2 A 2 X 2 +^4X4 (J%) +a 5 A 5 X 5 (l 5 ^\ ) = 0. 



f 5 (L)J 1 \f 5 (L)J 1 \f B (L) 

Here, f = (/l,/t), and f£ contain the electron distribution function parts 
in the reservoir k : as in Eqs. (jTUl fTT|) . Moreover, the matrices can be 
obtained using Eqs. (|T8|) — noting that the results for \E\ < A and \E\ > A 
are different due to the differing boundary conditions in the kinetic equations. 
Equation (|T9|) is linear, and the unknown fs(0) and fs(L) may be solved in a 
straightforward manner. Thus, for given temperatures and potentials in the 
reservoirs, we obtain the spectral currents ji and jx for each wire. 

To determine the thermopower, we have to find \i\ and [12 such that 
the condition ifT^! is satisfied. These potentials are expected to be small, so 
we can linearize the distribution functions in the reservoirs with respect to 
them, 



2k B T k J 2fc s T fe/ 

in the expressions for j Cj i and j Cj 2 obtained from Eq. lfT9)l . Thus, we obtain 
a linear equation for the potentials \x\ and //2, which can then be solved. 

However, without further approximations, the final expression for the 
thermopower is too complicated to provide much insight. We make the fol- 
lowing simplifications: 
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1. We approximate with its normal-state value T>i = 1. This is a good 
approximation only away from superconductor interfaces, but numerics 
show that it usually has little influence on the thermopower. 



2. Similarly, we approximate T>t ~ 1 in the latter term in Eq. (|18b|) . 
and retain it only in the first term. This is because Tand js provide 
more essential energy- dependent features than T>t, but we still wish to 
retain the temperature dependence of conductance due to the coeffi- 
cient in the first term. Comparing to exact numerics, we see that these 
approximations do not essentially affect the resulting thermopower. 

Using the above and the conservation of crAjs we can write Eqs. @ 
for the structure in Fig.^in a useful form at \E\ < A — noting that Aajs 
A^a^jsz = —A^a^iss according to the chosen directions: 



M" 



/ 


9i + 95 





-95 


^ 






-^ + gi (T 1 )+g 5 (T 5 ) 


5135 


^-95(T 5 ) 


-95 






-95 





92 +95 







V 


-^F-95(%) 


-95 


+ 92(T 2 ) + g 5 (T 5 ) 


5245/ 












(21a) 



M 



/7jvs(0)\ 

/t, 5 (0) 
h,5( L ) 

\fr,5(L)J 




(21b) 



where g\. = akA^/L^ is the conductance of wire k, = gi + gj + gk, 
9k = 9k/{^T^k) and {Tk) = Jo h % dx. For \E\ > A, the result is 



m: 



5135 

ft 
-95 



fll <T 1 > + <73<T 3 >+ 5 5<T5> ffi35 ^-g 5 (%) 



-55 \ 

Aajs 



2 



MI f = Cfi 2 + 



-95 

5245 

95(15) -55 92(^2) +9i(%)+g 5 (T 5 ) 5245/ 

(22a) 

/ 53 \ 

^+53<T 3 > 
54 

V -^s +g Jr 4 )J 



f° 

J L,3 
f 

J LA 



(22b) 



D 

Note that the temperature of the superconductors enters for energies \E\ > A 

,3/4- 



in the distribution functions /° . 
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Using Eq. lfT8)l with the assumption as above, we can calculate the cur- 
rents flowing in wires 1 and 2: 



(-AlGl3L,\\ 

-Aiaij T ,i 
\ A 2 a 2 jT,2 ) 



JC(MC-I)fi >2 , for \E\ < A 

|c(M A C-I)fi,2 + CM A Dfo, for|£|>A 



(23) 



Here, the sign of Jl/t,i 1S due to the choice of directions in Fig. ^ Now, 
we can integrate the rows 2 and 4 over the energy to obtain the observable 
charge currents, and linearize the distribution functions. Condition ({T^ then 
yields an expression of the form 



gdE 



L,l 
•0 

L,2 



J fJ +P3 (/jf J I 

(24) 

Here, g is a matrix containing temperature-dependent conductance-like quan- 
tities, and the integrals with yield, for example, terms proportional to the 
supercurrents flowing in the system at different temperatures. In the follow- 
ing, we simplify the expression further to obtain the dominant features. How- 
ever, Eq. l(24^l can well be used to obtain numerical thermopower estimates 
for given spectral data, without further approximations. 



L.O 

'0 

L.O 



3.6. Long-junction limit 

In this section, we discuss the dominant features of the thermopower, 
in particular the most significant contributions. We discuss first the long- 
junction limit A S> max[£r, ksT, //], and consider the effects of a finite A 
as corrections in the following sections. 

In the limit A 3> max[£ , r, ksT, fi] the induced potentials depend only on 
the temperatures T\ and T 2 of the normal reservoirs. This is due to the fact 
that the kinetic equations for \E\ < A are decoupled from the distribution 
function in the superconductors. In this limit, the spectral currents for \E\ > 
A are negligible, so they do not affect the integrated currents. This is because 
at these energies the distribution function is nearly constant in the system 
(Jl ~ 1, fr ~ 0), so there is no dissipative current. Additionally, at energies 
E ^> Et, the coefficient js ~ 0, so there is no non-dissipative supercurrent 
contribution in the current densities. 

If we take the limit A — > oo in expression l|24[l . approximate T« and 
neglect the energy dependence of T>t completely, we obtain the dominant 
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term in the induced potentials 

^ = ^(i? 1+ iCf^ 5 4 ) / ( 3 ^+l + i? 5 ) ^ 4 e{Is{Tl) ~ /S(T2)) ■ (25) 

Here, Is is the observable supercurrent Is{T) = Aa J^jsf^dE flowing in 
the system when the normal reservoirs are all at the same temperature T 
and there are no potential differences, and /° = tanh(i?/(2A;£T)). This 
expression states that a major contribution to the thermopower arises due 
to the temperature dependence 13 of the equilibrium supercurrent. 

Estimates for the equilibrium supercurrent appearing in Eq. (j23|l exist 
for SNS junctions. One is given in Ref. 33: 

eR SNS Is(T) = (k B T/E T fl 2 e -^V^T/E T gin(0) ^ (2g) 

where Et and Rsns = R3 + R4 + R5 are the Thouless energy and the 
resistance corresponding to the link between the superconductors. Expres- 
sion l[2T)|) is valid for Et,sns k B T <C A. Although the formula applies 
to a system where there are no extra normal-metal terminals, we can take 
them into account by scaling the supercurrent with the factor presented in 
Ref. 31: 

c i Asnsctsns ^ 

AsNS&SNS + 2(^1°"! + ^20"2) 

This holds for k B T S> Et,sns, Li,!® > Lsns, and £5 = 0. Obviously, in 
our case L5 / 0, but the results should still be useful, since the effect due to 
the central wire is small as discussed in Sect.I^Now, using Eqs. I)25|26|27jl . we 
obtain an estimate for the thermally induced potentials in the linear response 
limit Ti = T - AT/2, T 2 = T + AT/2: 



32 (A. - 2] vr 3 / 2 c scale ^3/4^5(2^4/3+^5) 

\V2 J scae (Ri + i?2 + R5)R 2 S ns (28) 
/27T - Z^E T /{k B T)) (k B T/E T ) e -^^/k B T/E T kBAT 



A*sc,l/2 



For example, in a structure where all the wires are similar, one would get for 
k B T = 3 E T 

Qns ~ 0.04^ « 3 ^— . (29) 
e K 

Comparing to numerical results (see Fig. O2J, estimate l(2*H)l is found to be 
useful for k B T > 3 Et. However, in the following we prefer to calculate the 
integrals with jg directly from numerical data, which gives better estimates 
at temperatures k B T < Et^sns- 
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Taking the coefficient Tinto account but still neglecting the energy de- 
pendence of T>t yields the following correction terms to the potentials: 

R f'°° 



% + R 5 )Rsns Jo 1 ' ' M 7 



T + i? 2 + fl 5 )#SiVS Jo 

Here, f^ k = tanh(E/ (2k bT^)). The correction is mostly small compared to 
the contribution from the supercurrent, but it is significant in some cases, for 
example at high temperatures UbT > 10 Et,sns- However, Tis not known 
analytically, so the integrals have to be calculated numerically. 

Thus, we have the following approximations for the potentials induced 
by a temperature difference: 

Hi ~ ^sc,l + ^T,l , ~ /Usc,2 + ^T,2 , (31a) 

Afj, = ii 2 — A*i ~ (^60,2 - Msc,i) + (mt,2 - Mr,i) • (31b) 

For left-right-symmetric structures, i?i = i?2, -^3 = ^4, these approximations 
yield /i^ / and Afi = 0, i.e. a finite N-S thermopower but a vanishing 
N-N thermopower. This is due to the symmetries T\{x) = — 72(^1,2 — x) and 
(75) = 0, which make the correction (fTTTTf) the same for both potentials n\ 
and [12 

However, there is also a third contribution, which arises from the energy 
dependence of T>t- This 7>r-effect has only a small impact on (Jns, but it is 
noticeable (see Fig. |4(b)| ) in Qnn in a left-right-symmetric structure. A fair 
estimate for the effect can be obtained by evaluating the matrix g in Eq. lffif|) 
using numerical results for T>t- For a left-right-symmetric structure, we can 
approximate the effect more roughly with 

r°°R 1 + 9 fsgfe / sech 2 ( ) sech 2 (^V)\ 

1 ^ 2R 3 +i? 5 v 7 

(32) 

where = {Vj^R^ are the energy-dependent spectral resistances of the 
wires. Generally this contribution is non-negligible only for large |T2 — T\\. 



3.7. Quasiequilibrium limit 

In the calculation above, we neglected the inelastic scattering. However, 
this assumption is not too restrictive for the thermopower, as long as a su- 
percurrent may flow through the normal metal. To see this, let us consider 
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a "quasiequilibrium" system, where inelastic scattering has relaxed the dis- 
tribution functions into Fermi functions with a local chemical potential n{x) 
and temperature T(x). Assuming the induced potentials are again small, ap- 
proximating T>t ~ 1, ~ 1, we can integrate the kinetic equations over the 
energy E (after multiplying the equation for by E). Using the fact that 
for quasiequilibrium distribution functions, in the limit where the induced 
potential differences are small, 



tix) = /"Mz) dE , T(x) = (j™E (1 - f L (x)) <1e) 2 , (33) 

we get the kinetic equations in wire i 

d x {Lid x yi{x) + eRiP s (T(x)) + Lif{x)d x T 2 ) = , (34a) 

d 2 x T 2 {x) = Q. (34b) 



Here 



r S (T(x)) = / f s (E)f L (E; T(x)) dE (35) 



is the supercurrent flowing in wire i. In the studied system, for a constant 
temperature, J| = — J| = /| = Is{T) and Ig = Ig = 0. Moreover, we may 
include the T^term through 

7W - \ JdET{E, r W ) _ J AE —^g-— , 

(36) 

but it does not have a direct physical interpretation. Both the spectrum 
js(E) of supercurrent-carrying states and the anomalous coefficient T(E) 
may in principle depend on the magnitude of inelastic scattering, but our aim 
here is to relate the measured thermopower to the temperature dependence 
of the actual supercurrent flowing in the system (and to the term T(x)). The 
supercurrent can be probed separately, and thus for this effect, we do not 
need to know the exact form of f s (E). 

The boundary and nodal conditions for the solutions fJ>i(x), Tj(x) are 
analogous to the general nonequilibrium case: in the reservoirs, they get the 
bulk values, except at the normal-superconducting boundary, say x = xns 
we have (as long as ksT <C A) 

d x T(x)\ x=XNS =0. (37) 



Moreover, at the nodes, the functions are continuous, and the charge and heat 
currents are conserved. These equations imply that the temperature T{x) is 
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constant in wires 3 and 4, and the heat current obeys the Wiedemann-Pranz 
law in the rest of the wires, i.e., T 2 (x) is a linear function of position. 

The solutions to these (linear) equations can easily be found. For sim- 
plicity, let us first ignore the term(s) 7\x). Denote the position in each wire 
by x, ranging from x = at the reservoirs or the left-hand node in wire 5, to 
Li at the other end. In wires 1 and 2, the solutions are 

fM = hi = const , fj, = fj,2 = const , (38a) 

T\x) = 2? + (T„ 2 3 - I?)f , T 2 (x) = T 2 + (T 2 4 - T|)f . (38b) 

Li L<2 

In wires 3 and 4, we get 

H(x) = hi— , //(x) = / u 2 — , (38c) 

T = T w3 , T = T W 4 , (38d) 



respectively. Finally, in wire 5 the solutions are 

eR 5 



x 



H(x) = Hi + eR 5 c 5 — + 

-^5 ^5 



I s (T(x'))dx', 



T\x) 



- w3 



ft 



wA 



l w3) 



L- 



(38e) 
(38f) 



Here [12, T-w3i T W 4, and C5 are constants that may be determined by 
requiring continuity and current conservation at the nodes. We get 



rp2 

1 w3 



(R 2 + R 5 )T 2 + RiTl 



R\ + R2 + R5 



rp2 

1 wA 



(Ri + R 5 )T 2 + R 2 T 2 



R-3 



C5 



Rsns 
Rsns 



R\ + R2 + R5 
RaIs{T w a) + (Ra + R 5 )Is(T w3 ) - R 5 



(39a) 



-4 I s {T)dT 2 



t's J wA 



T 2 

1 w3 



(39b) 



-{R 3 + R 5 )I s (T w4 ) + R 3 I s (T w3 ) + R 5 
Is(T w3 ), R S ns = R3 + Ra + R5- 



T ^ I s (T)dT 2 

Tl 3 T wA~ T w3 



(39c) 
(39d) 



Thus, the induced voltages Vf = Hi/ e ^i c = m the normal-metal 

reservoirs are determined from the observable supercurrent Is(T) between 
the superconductors, in the temperature range determined by 1^3 and T W 4. 
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Including the term T(x) is straightforward but leads to long expressions. 
Therefore, we present only the resulting voltages Vi\ 

Vi = V? + j^f- (T 2 4 - T 2 3 ) {%) + (T 2 3 - TfXTi) (40a) 

V2 = V? + (T 2 3 - T 2 4 ) {%) + (T 2 4 - T 2 )<T 2 ), (40b) 



Rsns 



where 



1 



Li 



[%) = - T^dx (41) 

L i JO 



is the average of the coefficient Tin wire i. 

Let us now take the linear-response limit around the temperature To, 
where 

AT 2 

T? /2 =T G 2 ±—, (42) 



and 

Tr,(T\ » Tr,(Tr,\ -L- 

dT 2 



I S (T)^I S (T ) + ^AT 2 . (43) 



In practice, the requirement for the validity of these formulae is AT 2 <C 
min(TQ , E^/kfj). In this case, always retaining only terms up to the first 
order in AT 2 , we get for the voltage in the left reservoir 

V sc = R 5 R 3 (2R i + R 5 ) dI s (T) ±T2 

1 2(R 3 + R i + R 5 )(R 1 +R 2 + R 5 ) dT 2 [ 1 

This is the same as the linearized form of Eq. (|2"K|) . Similarly, the linearized 
form of the Tcorrection follows that of Eq. l(30|) . This equivalence in the 
linear regime is due to the fact that the considered effect is not dependent on 
the exact shape of the distribution functions, and thus, in the linear regime, 
the relation between the thermopower and the temperature dependent equi- 
librium supercurrent is independent of the strength for inelastic scattering. 



4. GENERAL BEHAVIOR OF THE THERMOPOWER 

The typical magnitude and form of the temperature dependence of Qns 
is shown in Fig. [3 The figure shows that the magnitude is Qns ~ ^jr (often 
Qns,i ~ Qns,2)j and that the relevant energy scale is Et, corresponding 
to the link between the superconductors — which is also the energy scale 
of the spectral coefficients. Moreover, as shown in the preceding sections, 
the thermopower oscillates antisymmetrically with the phase difference (see 
Fig. OH. 
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5 10 15 



k B T l /E T 

Figure 2. Temperature dependence of the N-S thermopower in a left-right- 
symmetric structure — Lj,&j and Aj assumed equal. Phase 4> = tt/2 is 
fixed. 



The N-S thermopower shows re-entrant behavior: it vanishes as Xi, — > 
0, peaks at an intermediate temperature, and decays as T\,T<i — ► oo. The 
peak near T ~ 2Ex/kB corresponds to the position where Is{T) changes 
rapidly as a function of the temperature, as Eqs. (|25J26|) indicate. Moreover, 
Eqs. Il2. r )|:-i()[i predict that Qns decays at least as T~ 2 for T\ w T2. In fact, 
Eq. l{28|) implies that the contribution from the supercurrent should decay 
exponentially, but the contribution from Thas a slower speed of decay. 

The most significant contribution to Qns is described by Eq. (|25[l: it 
yields almost all of the effect seen in Fig. [21 The correction due to Taccounts 
for the most of the deviation from the numerical result, and has an effect be- 
low 10% in relative magnitude for temperatures T\,T2 < 10 Er/kB, but the 



proportion becomes relatively larger at higher temperatures — see Fig. 5(a) 



Moreover, T>t has a negligible effect on Qns, as seen in Fig. 5(a) 



The thermopower Qns with respect to the actual (non-equilibrium) su- 
percurrent flowing in the structure is shown in Fig. 03 For a given tempera- 
ture, the relation between Qns an d Is appears to be linear, in accord with 
Eq. (EE). 

The temperature dependence of Qnn shown in Fig. |4(b)| is very dif- 
ferent from that of Qns, since the structure there is assumed to be left- 
right symmetric, so that for Qnn the contributions from Is and T vanish 
and the Pr-effect dominates (see Fig. |5(bjj ), In this case, Eq. l(32")l predicts 
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Figure 3. Oscillation of the thermopower. Left: Phase-oscillation of the 
thermopower at T2 « T\ = 1 Et- Right: Relation between Qns and the 
non-equilibrium supercurrent Is at different temperatures T = T\ « Ti- 
Magnitude of the supercurrent here is nearly the same as in equilibrium. 
Due to a numerical convergence problem, data for phases |0| > 2.1 is not 
available, and hence substituted with straight line segments (dotted). In 
the right picture, the dotted part also approximately describes the accessible 
non-dissipative regime in the case when the supercurrent is externally driven. 




(a) (b) 

Figure 4. Temperature dependence of the N-N thermopower. Left: Qnn in 
an asymmetric structure, R3/2 = R\ = R2 = R4 = R5 with Ao equal, at 
<fi = 7r/2. Right: Qnn in a left-right-symmetric structure, where the T>t- 
contribution dominates. The thick white lines indicate a change of sign for 
R{T\) — R(T2), which occurs near the change of sign of Qnn- 
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(a) Qns,2 at Ti « T2 



(b) Qnn with fixed Ti = 3.6 Et 



Figure 5. Thermopower at (f> = tt/2 in a symmetric structure, L and Ao 
assumed equal for each wire. The numerical results and different analytic 
approximations from Sect. IH.fil are shown. Here, "Is (analytic)" is the ap- 
proximation from Eq. 1281 



Afi oc (Ti - T 2 ) 2 at T\ ~ T2 in symmetric structures, thus the thermopower 
appears only at large temperature differences, i.e. in the non-linear regime. 
For asymmetric structures, 2>r~effect is less important, as Eq. (|31b|) implies 
A/i oc (Ti — T2) at T\ T2, and this effect tends generally to wash out the 
I?T-effect also at large temperature differences. This behavior can be seen 



in Fig. 4(a) where Qnn has a similar temperature dependence as Qns in 
Fig. 10) ■ Effects due to the sample geometry are discussed in more detail in 
the following section. 

Note that the re-entrant temperature-dependent behavior of the con- 
ductance 14 induced by the proximity effect also arises due to the energy de- 
pendence oiT>T- Thus, the finite NN-thermopower in a left-right-symmetric 
structure may at least qualitatively be understood to be caused by an in- 
duced left-right asymmetry in the resistances, which makes the potential 
induced in the colder electrode different from the one in the hotter electrode. 
Indeed, the sign changes in Fig. |4(bj| (along the diagonal and the 1/Ti-like 
curve) occur close to the curves where the (zero-bias) conductances, 16 

c|T) .fr.?« dE . jj = <Bf . >fl , (45) 







of wires 1 and 2 at the temperatures T\,Ti become equal. Approximately 
such behavior is also expected on the basis of Eq. (pS2"|h 
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(a) Qns,2 (b) Qnn 

Figure 6. Thermopower in an asymmetric structure with all except one wire 
of length L, at T\ ~ T2, (ft = ir/2. Aa assumed equal in all wires. 



Comparing to numerical results, we find that the approximations ob- 
tained from Eq. (|2"3|) in the long-junction limit are quantitatively accurate 
in a symmetric system, both for Qns an d Qnn (the relative error in Qns 
is less than 5% nearly everywhere in Fig. although near T\,T2 ~ it is 
~ 15%, at worst). In setups geometrically deformed from this, the accuracy 
may decrease slightly, but qualitative features are still well retained (in all 
cases tested). Equations ()24l2. r )|:-i()|:-i2ll may thus be used to provide use- 
ful approximations to the induced potentials, if the spectral coefficients are 
known. 



5. DEPENDENCE ON SAMPLE GEOMETRY 

The thermopower induced by the proximity effect has a non-trivial de- 
pendence on the exact geometry of the structure: varying the geometry of 
the structure changes both the magnitudes and the temperature dependence 
of the induced potentials, as seen in Fig. The changes arise partly in the 
kinetic equations, the effect of which is estimated by the pref actors in the 



approximations l|25l30p . but also the spectral equations contribute through 
the geometry dependence of the spectral coefficients — their behavior is dis- 
cussed in Sect. E] 



An effect visible in Fig. 6(a) is that compared to the symmetric structure, 



the induced potentials typically tend to be smaller in structures where the 
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lengths of some of the wires are strongly different from each other. This be- 
havior appears due to both the spectral coefficients and the kinetic equations. 
For example letting -^1,2/^3,4,5 — ► tends to decrease the thermopower, as 
then Is and (T) are suppressed, while in the opposite limit most of the tem- 
perature drop occurs in wires 1 and 2 where there is little coupling between 
the charge and energy currents, leading again to a smaller thermopower (as 
indicated by the pref actors in Eqs. I)25l30|l ). 

The effect the geometry has on the induced potentials is more significant 
for Qnn than for Qns, since the contributions to A/i depend strongly on the 
asymmetry and may also vanish due to symmetries in the structure. Gen- 
erally, the approximations imply that in left-right-asymmetric structures, 
Qnn should be finite and roughly resemble the N-S thermopower in temper- 



ature dependence, as in Fig. 4(a) Moreover, this thermopower should be 
discernible even for small amounts (~ 10%) of asymmetry in resistances, as 
can be seen in Fig. |6(b)[ 

Although the temperature dependence of thermopower generally scales 
with Et,sns similarly to the spectral coefficients, its exact form depends on 
the geometry of the structure. That is, the terms due to js and (T) appearing 
in the thermopower estimates vary in the geometry dependence, hence one 
can alter their relative weights by adjusting the structure. Moreover, they 
may also differ in sign and in the way they depend on the temperature, 
so their weighing has a direct effect on the temperature dependence of the 
thermopower. This is important especially for the N-N thermopower, which 
can be seen in Fig. |6(b)| where changes in L\ and L3 cause qualitatively 
different results. 



5.1. Thermopower and the central wire 

As found analytically in Sect. I.3.6.( the spectral coefficients that affect 
the thermopower most are js and T. In the long-junction limit, js couples 
the energy and charge currents in the central wire, and (^1/2) do this in the 
wires leading to the N-reservoirs. By changing the proportion of the lengths 
of these wires, using the factors in Eqs. (|25|H()jl as a guide, the weaker effect 
due to (T) may be brought to dominate. In the following, we consider a 
special case where T should manifest: L5 <C Lsns- (Another possibility 
suggested by Eqs. (j2SI3flj) would be L\,L 2 S> Lsns, but in this case, the 
terms (^1/2) become small.) 

In the limit L5 / min[Li, L 2 , L3, L4] — ► 0, the supercurrent-effect approx- 
imation ({25} yields no contribution, and the second term in Eq. tpTTTfl vanishes, 
leaving only the contribution due to (7^/ 2 )- If the structure is left-right- 
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Figure 7. N-S thermopower for a structure with the central wire short or 
absent, L\ = L 2 = L4 and Aa assumed equal. Moreover, T\ « T 2 and 
<j> = 7r/2. The "full analytic approximation" is obtained using Eq. (|24"|) and 
numerical data for T>t, js and T. The approximation for the contribution 
due to Tobtained from Eq. (dashed) is nearly the same in absolute value 
for all cases shown. 



symmetric, then T\/ 2 = 0, and the thermopower also vanishes (confirmed by 
simulations). In fact, this can be shown exactly by noting that fx = is the 
exact solution to the kinetic equations in this special case. 

In asymmetric structures especially with L3 ^ L4, generally (^1/2) 7^ 
and their sign is the same for both wires 1 and 2. Hence Qns becomes 
finite and there is a large difference between Qns,i and Qns,2, as the con- 
tributions from Eq. P?)]) differ in sign for electrodes 1 and 2. In Fig. the 
approximation lj30|l matches the numerical result well: and (T2) yield 

the dominant contribution. 

As the terms l(25|l due to Is have the same sign for both electrodes and 
vanish as L5 — > 0, Qns may change its sign with increasing temperature for 
intermediate length L5 of the central wire. This is possible because for a 
certain range of values for L5, Is contributes the same overall amount as 
(Ti/ 2 ), but in the opposite direction for either Qns,i or Qns,2- As the two 
contributions have slightly different temperature dependence, this can result 
in a change of sign in the thermopower, for a suitable L5. The effect is 
illustrated in Fig. d where a change of sign occurs for L5 = 0.1L. This sign 
change may partially explain the results in Ref. 21, see Subs. IfTTI 
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(a) Qns, To ~ Ti ~ T 2 (b) Qnn, To = Ti = 3.6 St 

Figure 8. Thermopower for different A/Et- Left: Qns i n a structure with 
all wires equal. Right: Qnn in a structure with L\ = L2 = L4 = L5 = L and 
L3 = L or L3 = 2L. In both figures, Aa is assumed equal for all wires. The 
temperature dependence of A is neglected in the figures. 



5.2. Effect of decreasing A/E T 

In reality, the superconducting energy gap A may be quite small in 
systems realized experimentally, and not necessarily much greater than Et- 
For example, Eq. (fT^jl shows that A = 2 K corresponds roughly to 10 Et 
for Lsns = 1 A* m an d to 3 Et for Lsns = 0.5 /mi, for a typical value 
D = 200 ssi of the diffusion constant. 

For a finite A, the spectral currents from the energies E > A contribute 
to the integrated currents. This in turn couples the temperature To of the 
superconductors to the system, as seen in Eq. (|24|l. Additionally, the spectral 
coefficients at energies E < A are modified from their long-junction limits, 
and hence the contribution from E < A is also affected. 

Estimates for the effect due to the coupling of Tq may be obtained from 
Eq. Generally, terms proportional to (/° — /° {), — /° 2 ) and 

(/l 1 — appear for E > A, thus any of the corresponding tempera- 

ture differences induces potentials. However, these contributions tend to be 
smaller than those from E < A, at least for A/Et > 10 — supposing of 
course that T\ — T2 is of the order T\ — Tq or T2 — Tq, whichever is larger. 



Figure 8(a) shows the effect of decreasing A/Et on the N-S thermopower. 



The energy scale of the temperature dependence seems to be compressed to- 
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wards T = 0, but no significant qualitative changes occur. The result is 
similar if T\ = Tq are kept fixed and T2 is varied. 

As Qnn is generally smaller than Qns, the relative changes tend to be 
larger in it. This can be seen in Fig. |8(b)| where the contribution from E > A 
for A/Et = 10 peaks near T2 = 6 in the left-right-symmetric structure, 
and is also the dominant effect there. For the asymmetric structure, the 
change is relatively much less significant. 

6. COMPARISON WITH RECENT WORK 

In this section, we compare the theory presented in this paper and in 
Ref. 26 to recent experiments and other theories on the thermopower in 
proximity effect structures. 

6.1. Comparison to experiments 

To our knowledge, there exist so far five published experiments from 
two groups probing the thermopower in Andreev interferometers 18-22 . One 
can distinguish three features of these experiments to which our theory can 
be compared: magnitude of the observed effect and the dependence on the 
temperature and magnetic flux. 

In Ref. 18, the temperature gradient produced by the heating current 
was not explicitly measured, but only estimated using a heating model. Such 
a procedure yielded too large an estimate (of the order of //V/K) for the 
measured NN-thermopower. 20 This observation is supported by Refs. 19,20, 
which discuss further measurements by the same group, now with the direct 
measurement of the temperature. Below, we concentrate on the results of 
Ref. 19 — the observations in Ref. 20 are quite analogous. These measure- 
ments yield a thermopower of the order of 50-60 nV/K at the temperature of 
the order of one Thouless temperature corresponding to the distance between 
the two superconductors. The measured structure slightly differed from that 
considered in our manuscript, as the superconductors were fabricated directly 
on top of the normal wires. However, we may assume the N-S contact resis- 
tances to play a similar role as the resistances of the wires 3 and 4 (see also 
Subs. lo~2~|) , We also note that the magnitude of the NN thermopower greatly 
depends on the asymmetry in the structure. For example, taking the follow- 
ing estimates for the resistances of the five wires: R\ = 5 O, R2 = 5.5 f2, 
R3 = Ri = 0.5 f2, i?5 = 5 O (note that mostly only the ratios of the re- 
sistances are relevant) and assuming Ex/ks = 275 mK (corresponding to a 
wire of length 700 nm) we get Q NN = 60 nV/K at T = 295 mK, 4> = it/2. 
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This is very close to the experimental results. 

Concerning a comparison to the experiments in Ref. 21, we obtain NS 
thermopower which is larger by an order of magnitude (of the order of 2 
/iV/K, whereas the experiments report some 70 nV/K). However, it is not 
totally clear to us where the normal-metal reservoirs in this experiment re- 
side, and hence, which values should be used for the resistances R\ and i?2- 

Our theory predicts a nonmonotonic thermopower as a function of the 
lattice temperature, with a maximum at around Et/Isb, given by the Thou- 
less energy corresponding to the wire length between the superconductors. 
This is in agreement with Fig. 4 in Ref. 18. The temperature scale of the 
dependence Q(T) in Refs. 21,22 is of the order of Et/^b, but these pa- 
pers also report a sign reversal of the thermopower as a function of T. As 
indicated in Subs. I5.1.( this may be possible in a suitable (left-right asym- 
metric) geometry where the overall supercurrent and 7^term contributions 
to the thermopower are comparable. Another possible source for the sign 
change is the additional term induced by the coupling of the superconductor 
temperature Tq at high temperatures. 

As shown in Subs. 13.4.1 the exact solution of our equations (within the 
assumption of electron-hole symmetry) is an antisymmetric function of the 
flux (phase) with respect to zero flux. This is in agreement with the sym- 
metries obtained for the "Parallelogram" structure in Ref. 18, and those 
measured in Refs. 19-22. However, our results cannot explain the symmetric 
flux dependence seen in the "House" interferometer in Ref. 18. 



6.2. Comparison with other theories 

The first theoretical discussion of thermoelectric phenomena under the 
superconducting proximity effect by Claughton and Lambert 17 showed how 
the thermoelectric coefficients can be calculated from the scattering theory. 
The general qualitative features one can extract from the analytic formulae 
are valid in any kind of structures composed of normal metals and super- 
conductors. However, the numerical simulations of these systems (included, 
for example, in Refs. 17,23) suffer from the small size of the simulated struc- 
tures which makes it difficult to differentiate between the effects related with 
electron-hole asymmetry (large in the small simulated structures) from the 
dominant effects in realistic experimental samples with thousands of chan- 
nels. 

In Refs. 24,25, the thermopower of normal-superconducting structures 
is considered in the limit of a weak proximity effect (with large NS interface 
resistance) and in the linear regime. To compare our theory to those in their 
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Figure 9. The spectral supercurrent js for different ratios of wire lengths. 



f £ * cosh "(dk) iE - (46) 



works, let us set T 2 = T + 5T/2, T x = T - ST/2 in Eq. for the S-N 
thermopower in a left-right symmetric structure, 

Q = ^2 = 1 ^5 1_ 

y S N - §T - 4 2jRi ^ + ^ (A;bTo)2 

Here, j' s = L^js- This is analogous to Eq. (7) in Ref. 24 and to the upper 
line of Eq. (8) in Ref. 25. The primary coefficient g z+ in these equations is 
proportional to the spectral supercurrent across the normal-superconducting 
interface, analogous to f s in Eq. (|4*6|) . 

The lower line of Eq. (8) in Ref. 25 describes a thermopower arising 
between the normal-metal reservoirs. This term arises from the spectral 
supercurrent across the interface at the energies above the superconducting 
gap A and is thus similar to the contributions discussed in Subs. 15.2.1 As 
shown there, these contributions are mostly important only for highly left- 
right symmetric structures in the case when L$ns is n °t much larger than 
the superconducting coherence length, or when T\ or T 2 are close to A/ks- 



7. BEHAVIOR OF SPECTRAL COEFFICIENTS 

Since the thermopower is induced by the spectral coefficients js and (7), 
we describe here briefly their dependence on the energy and on the sample 
geometry. 



7.1. Spectral supercurrent 



The spectral supercurrent js is the better-known of the two main con- 
tributors for the thermopower. In structures such as considered here, it typ- 
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Figure 10. Behavior of Tin a structure with L\ = L 2 = L$/l.l = L 4 = L 5 
and Aa equal, at <\> = n/2. Left: The energy-dependence of the averaged 
dimensionless coefficient (T) = J^ k Tdx in different wires. Right: Spatial 
behavior of J^°TdE / Et- Different branches correspond to T in different 
wires. Locations of the two nodes are indicated by the dotted lines, and 
reservoir interfaces are denoted by dots. 



ically has the energy dependence shown in Fig. |HJ the characteristic energy 
scale is Et,sns, and the dimensionless quantity LsnsJS does not depend 
greatly on the exact geometry of the structure (supposing L\,L2 > Lsns)- 
The behavior of js is discussed in detail in Ref. 31, and the results are also 
mostly applicable to the structure considered here, as the central wire typi- 
cally causes no significant qualitative changes. 



7.2. Anomalous T-term 

The T-term l|8d|l . whose effect is neglected in most previous papers, is 
a coefficient appearing in the non-equilibrium part of the supercurrent, and 
it has no clear physical interpretation. Nonetheless, it affects the potentials 
induced due to a temperature difference, and provides in some cases the 
dominant contribution which may be large by itself. Hence it is interesting 
to examine the energy and spatial behavior of this coefficient, which is here 
considered in the long-junction limit A 3> Et- 

The typical energy dependence of the average Tin different wires of the 
structure in Fig. ^ is shown in Fig. Its magnitude is approximately less 
than 10% of the spectral supercurrent LsnsJS, but the energy dependence is 
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rather similar: a single peak appears on the scale of E^^ns although there 
is less oscillation than in LsnsJS- 

The spatial dependence shown in Fig.ESgives insight to the magnitudes 
of Tin different wires. Here, Tvanishes at the reservoirs and changes its sign 
near the center of the SNS link, forming peaks between the S-reservoirs and 
the center. This leads to (T5) in nearly left-right-symmetric structures, 
as seen in Fig. Moreover, T3 and T are typically the largest and the two 
differ in sign. 

In the wires 1 and 2 leading to the N-reservoirs, T decays monotonically 
to zero on the length scale of Lsns (according to numerical results). In 
fact, no sign changes are possible for 75y 2 , as x = const(x) in these wires. 
Moreover, note that for the quantities (T/2) relevant for the thermopower, 
this kind of a decay would imply that (T/2) vanish for Li^/ Lsns ~ ► 00 • 

The magnitude and the sign of T in the wires 1 and 2 depend on the 
points where the wires are connected to the SNS link, since Tis continuous 
in the structure. Moreover, if the N-reservoirs are far {L\,L2 > Lsns)-, the 
choice of the connection points does not affect much the Tin the SNS link, 
hence one may estimate the qualitative behavior of T in wires 1 and 2 by 
looking at Fig. EB 

Using the symmetry arguments from Subs. I3.4.| one finds that T os- 
cillates antisymmetrically with the phase difference. Moreover, the overall 
phase behavior of Tseems to be similar to that of js, according to numerical 
results. 



8. DISCUSSION 

In this paper, we have shown how the presence of the supercurrent can 
lead to a finite thermopower in an Andreev interferometer even in the pres- 
ence of complete electron-hole symmetry. Generally, the thermopower probes 
the temperature dependence of the supercurrent, which in the long-junction 
limit is determined by the Thouless energy corresponding to the distance be- 
tween the superconductors. The actual magnitude of the thermopower also 
strongly depends on the studied geometry of the sample, which has to be 
taken into account when comparing to the experiments. 

Another interesting result is the thermopower induced by the "anoma- 
lous" coefficient T This coefficient has been identified in a number of the- 
oretical publications (for example, Refs. 24,25,27,34-36 — it is frequently 
referred to as the "anomalous current" as its properties are not very well 
known) but typically its effect is neglected and we are not aware of any other 
observable whose behavior would be dictated by T 
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The theory presented in this paper and in Ref. 26 seems to be in fair 
agreement with most of the experimental results published so far (see Subs. 
Ifi.l.jl . However, unless device geometry (resistances of the wires) are quite 
well known, it is difficult to make a detailed comparison of the magnitude 
and temperature dependence of the measured thermopower. There exists 
also one experimental result (the symmetric oscillations of the thermopower 
in the "House" interferometer in Ref. 18) which clearly cannot be explained 
with the present theory. 

Confirming the theory clearly calls for more experiments on the sub- 
ject. Rather than as a function of the applied flux, one could probe the 
thermopower by externally applying a supercurrent between the supercon- 
ductors (c.f., Fig- EJ - In this way, one would also be able to compare the ob- 
served thermopower to the temperature dependence of the equilibrium (crit- 
ical) supercurrent, essentially with no fitting parameters — the normal-state 
resistances can be fairly accurately measured in the considered multiprobe 
structure. Further, it would be interesting to confirm the relation between 
the thermopower and the "anomalous" term T, either by substracting the su- 
percurrent part from the results or, for example, in a structure without the 
central wire. The effect of the additional terms, present in short wires (due to 
a finite A/Ej<) or highly symmetric structures (due to the proximity-induced 
temperature dependence in the resistances), would also be worth studying, 
but perhaps more difficult to realize in practice. 



ACKNOWLEDGMENTS 

We thank Frank Hekking, Nikolai Kopnin, Mikko Paalanen and Jukka 
Pekola for useful discussions. 



REFERENCES 

1. V. L. Ginzburg. Nobel lecture (2003). 

http: / /www. nobel.se /physics /laureates / 2003 /ginzburg- lecture. html 

2. V. L. Ginzburg. Zh. Eksp. Teor. Fiz., 14, 177 (1944). [J. Phys. USSR 8, 148 
(1944)]. 

3. N. W. Ashcroft and N. D. Mermin. Solid-state physics. Saunders College Pub- 
lishing (1967). Chapters 12-13. 

4. D. J. van Harlingen. Physica, 109 & HOB, 1710 (1982). 

5. Y. M. Galperin, V. L. Gurevich, V. I. Kozub, and A. L. Shelankov. Phys. Rev. 
B, 65, 064531 (2002). 

6. C. J. Pethick and H. Smith. Phys. Rev. Lett., 43, 640 (1979). 

7. J. Clarke, B. R. Fjordb0ge, and P. E. Lindelof. Phys. Rev. Lett., 43, 642 (1979). 



Thermopower in Andreev Interferometers 



8. A. Schmid and G. Schon. Phys. Rev. Lett, 43, 793 (1979). 

9. J. Clarke and M. Tinkham. Phys. Rev. Lett., 44, 106 (1980). 

10. C. J. Lambert and R. Raimondi. J. Phys.: Condens. Matter, 10, 901 (1998). 

11. H. Courtois and B. Pannetier. J. Low Temp. Phys., 118, 599 (2000). 

12. W. Belzig, F. K. Wilhelm, C. Bruder, G. Schon, and A. D. Zaikin. Superlatt. 
Microstruct., 25, 1251 (1999). 

13. P. Dubos, H. Courtois, B. Pannetier, F. K. Wilhelm, A. D. Zaikin, and G. Schon. 
Phys. Rev. B, 63, 064502 (2001). 

14. P. Charlat, H. Courtois, Ph. Gandit, D. Mailly, A. F. Volkov, and B. Pannetier. 
Phys. Rev. Lett., 77, 4950 (1996). 

15. Yu. V. Nazarov and T. H. Stoof. Phys. Rev. Lett., 76, 823 (1996). 

16. E. V. Bezuglyi and V. Vinokur. Phys. Rev. Lett., 91, 137002 (2003). 

17. N. R. Claughton and C. J. Lambert. Phys. Rev. B, 53, 6605 (1996). 

18. J. Eom, C.-J. Chien, and V. Chandrasekhar. Phys. Rev. Lett., 81, 437 (1998). 

19. D. A. Dikin, S. Jung, and V. Chandrasekhar. Phys. Rev. B, 65, 012511 (2002). 

20. D. A. Dikin, S. Jung, and V. Chandrasekhar. Europhys. Lett, 57, 564 (2002). 

21. A. Parsons, I. A. Sosnin, and V. T. Petrashov. Phys. Rev. B, 67, 140502 (2003). 

22. A. Parsons, I. A. Sosnin, and V. T. Petrashov. Physica E, 18, 316 (2003). 

23. T. T. Heikkila, M. P. Stenberg, M. M. Salomaa, and C. J. Lambert. Physica B, 
284-8, 1862 (2000). 

24. R. Seviour and A. F. Volkov. Phys. Rev. B, 62, 6116 (2000). 

25. V. R. Kogan, V. V. Pavlovskii, and A. F. Volkov. Europhys. Lett, 59, 875-881 

(2002) . 

26. P. Virtanen and T. T. Heikkila. Phys. Rev. Lett. 92, 177004 (2004). 

27. T. T. Heikkila, T. Vanska, and F. K. Wilhelm. Phys. Rev. B, 67, 100502(R) 

(2003) . 

28. Y. V. Nazarov. Superlatt. Microstruct, 25, 1221 (1999). 

29. J. Rammer and H. Smith. Rev. Mod. Phys., 58, 323 (1986). 

30. K. D. Usadel. Phys. Rev. Lett, 25, 507 (1970). 

31. T. T. Heikkila, J. Sarkka, and F. K. Wilhelm. Phys. Rev. B, 66, 184513 (2002). 

32. W. Hansen. Phys. Rev. B, 31, 3504 (1984). 

33. F. K. Wilhelm, G. Schon, and A. D. Zaikin. Phys. Rev. Lett, 81, 1682 (1998). 

34. S.-K. Yip. Phys. Rev. B, 58, 5803 (1998). 

35. S.-W. Lee, A. V. Galaktionov, and C.-M. Ryu. J. Korean Phys. Soc, 34, S193 
(1999). 

36. E. V. Bezuglyi, V. S. Shumeiko, and G. Wendin. Phys. Rev. B, 68, 134506 
(2003). 



